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<N ; Abstract 

The magnetic moment of the Roper resonance is calculated in the framework of a low-energy 

CN ■ effective field theory of the strong interactions. A systematic power-counting procedure is imple- 

t>X) mented by applying the complex-mass scheme. 
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I. INTRODUCTION 



Chiral perturbation theory [EH] provides a successful description of the Goldstone boson 
sector of QCD (see, e.g., Ref. |3| for a recent review). A straightforward power counting, i.e. 
correspondence between the loop expansion and the chiral expansion in terms of momenta 
and quark masses at a fixed ratio [2], is obtained by using dimensional regularization in 
combination with the modified minimal subtraction scheme. Therefore, a systematic and 
controllable improvement is possible in perturbative calculations of physical quantities at 
low energies. The construction of a consistent power counting in effective field theories with 
heavy degrees of freedom turns out to be a more complex problem. For example, power 
counting is violated in baryon chiral perturbation theory if dimensional regularization and 
the minimal subtraction scheme are applied jij]. The problem has been handled by employing 
the heavy-baryon approach [5| and, alternatively, by choosing a suitable renormalization 
scheme p-lSl]. Using the mass difference between the nucleon and the A(1232) resonance 
as an additional expansion parameter, the A resonance can also be consistently included in 
the framework of effective field theory [l0-14|. On the other hand, the inclusion of heavier 



baryon resonances such as the Roper resonance requires a non-trivial generalization. In 
this case t he p roblem of power counting can be solved by using the complex-mass scheme 
(CMS) [l5-19 which can be understood as an extension of the on- mass-shell renormalization 
scheme to unstable particles. In previous papers we have calculated the pole masses and 



the widths of the p meson and the Roper resonance |20l . |21| . In the current paper we 



consider the magnetic moment of the Roper up to 0(q z )} While the extraction of these 
quantities from experimental measurements at present seems to be unfeasible, our expression 
for the magnetic moment may be used in the context of lattice QCD. Effective field theories 
predict the quark-mass dependence of physical observables and can be used to extrapolate 
simulations in the framework of lattice QCD performed at unphysically large masses of the 
light quarks. In return, lattice QCD provides a way to determine the low-energy constants 
from the underlying theory. 



II. EFFECTIVE LAGRANGIAN 



In this section we specify the effective Lagrangian relevant for the subsequent calculation 
of the electromagnetic vertex of the Roper at 0(q 3 ). We include the pion, the nucleon, the 
Roper, and the A as explicit degrees of freedom. The effects of other degrees of freedom are 
buried in low-energy coupling constants. We write the effective Lagrangian as 2 



C — Co + Ar + Cr + £nr + £ 



AR ■ 



(1) 



where Cq is given by 

C = N(iJp-m m )N + R(ilp-m R o)R 



(ilp - m A0 ) g 1 " - % (YD U + 1 U D^) + % Yfl Y + ™ao 



F*u- (2) 



1 Here, q stands for small parameters of the theory such as the pion mass. 

2 To simplify the notation only bare masses are supplied with a subscript 0. 
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Here, N and R denote nucleon and Roper isospin doublets with bare masses m^o and ?bro, 
respectively. ^ v represents the vector-spinor isovector-isospinor Rarita-Schwinger field of 
the A resonance 22j with bare mass mAo, is the isospin-3/2 projector (see Ref. 14J for 
more details). The covariant derivatives are defined as follows: 



r„ 



1 

2 



(3) 



where H stands either for the nucleon or the Roper. The pion fields are contained in the 
unimodular, unitary, (2x2) matrix U and u = y/U. The external electromagnetic four- vector 
potential A^ enters into t> M = — e ^ A^ and vfi = — § A^ (e 2 / (An) 1/137, e > 0). 

The lowest-order Goldstone-boson Lagrangian including the quark-mass term and the 
interaction with the external electromagnetic four-vector potential A^ reads 



F 2 F 2 M 2 F 2 

Hf = Tr (dJJdPXP) + -^—T? + U)+i — Tr 



(4) 



F denotes the pion-decay constant in the chiral limit: F n = F[l + 0(q 2 )] = 92.4 MeV; M is 
the pion mass at leading order in the quark-mass expansion: M 2 = 2Brh, where B is related 
to the quark condensate (qq)o in the chiral limit 0. 

The interaction terms Cr, Cnr, and £a_r are constructed in analogy to Ref. [23j. The 
leading-order (0(q)) pion-Roper coupling is given by 



C 



where gn is an unknown coupling constant and 

u^ — i u^d^u — ud^iv — i 
The second- and third-order Roper Lagrangians relevant for our calculation read 



(5) 



(6) 



where 
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(3) 
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+ H.c. +2id* 7 R £>^R + H.c. + 



(7) 
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and Cg, Cy, tig, and <iy are unknown coupling constants. The ellipsis denote those terms of 
the most general second- and third-order Roper Lagrangians which do not contribute to the 
electromagnetic vertex of the Roper at 0(q 3 ) and H.c. refers to the Hermitian conjugate. 
The leading-order interaction between the nucleon and the Roper is given by 



(i) 

NR 



9NR 



RYlsu^N + H.c. 



(9) 
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with an unknown coupling constant gNR- Finally, the leading-order interaction between the 
A and the Roper reads 

4* = -9ar ^ d (<T + z YY) u v R+ H.c. , (10) 

where g^R is a coupling constant and we take the "off- mass-shell parameter" z = — 1. Note 
that at 0(q 3 ) the NA La grangian does not contribute. 



III. PERTURBATION THEORY, RENORMALIZATION, AND POWER COUNT- 
ING 



The CMS |15Hl9l] originates from the Standard Model where it was developed to derive 
properties of W, Z , and Higgs bosons obtained from resonant processes. What makes the 
situation somewhat different in the case of the strong interactions is the fact that hadrons, 
including resonances, are thought to be composite objects made of quarks and gluons. The 
characteristic properties of hadron resonances eventually have to be described by QCD. 
Within the present effective-field-theory approach, to a given resonance we assign an ex- 
plicit field with corresponding spin, isospin, and parity content. Furthermore, for a generic 
resonance R, we introduce a complex renormalized mass z x defined as the location of the 
corresponding complex pole position in the chiral limit, z x = rriR X — iTR X /2. We assume 
Tr x to be small in comparison to both vn^ x and the scale of spontaneous chiral symmetry 
breaking, A x = AttF. Corrections to the complex pole position due to the finite quark 
masses are treated perturbatively. Our perturbative approach to EFT is based on the path 
integral formalism. In this framework the physical quantities are obtained from Green's 
functions represented by functional integrals. The integration over classical fields corre- 
sponding to particles is performed in the standard way, i.e., the Gaussian part is treated 
non-perturbatively and the rest perturbatively. In particular, the functional integral is per- 
formed for both stable and unstable degrees of freedom. For stable particles the path integral 
formalism is equivalent to the operator formalism based on the Dirac interaction representa- 
tion. Unfortunately, it is not obvious how to apply this representation to field operators for 
unstable particles, because, strictly speaking, there is no free Hamiltonian for unstable par- 
ticles. Therefore, we stick to the functional integral where one can perform the integration 
independently of the nature of the field. 

In the following, we apply the CMS to have a consistent power counting also applicable 
to loop diagrams. This renormalization scheme is realized by splitting the bare parameters 
(and fields) of the Lagrangian into, in general, complex renormalized parameters and counter 
terms. We choose the renormalized masses as the poles of the dressed propagators in the 
chiral limit: 

m R0 = z x + Sz x , 
m m = m x + 5m , 

m A o = za x + $ z a x , (11) 

where z x is the complex pole of the Roper propagator in the chiral limit, m x is the mass of 
the nucleon in the chiral limit, and z& x is the complex pole of the A propagator in the chiral 
limit. We include the renormalized parameters z x , m, and z& x in the propagators and treat 

the counter terms perturbatively. The renormalized couplings Cg and of are chosen 
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such that the corresponding counter terms exactly cancel the power-counting-violating parts 
of the loop diagrams. 

While the starting point is a Hermitian Lagrangian in terms of bare parameters and 
fields, the CMS involves complex parameters in the basic Lagrangian and complex counter 
terms. Although the application of the CMS seems to violate unitarity, the bare Lagrangian 
is unchanged and unitarity cannot be violated in the complete theory. On the other hand, it 
is not obvious that the approximate expressions to the S'-matrix generated by perturbation 
theory also satisfy the unitarity condition since the conventional Cutkosky cutting equations 
(241 are not valid in the framework of CMS. However, it is possible to derive generalized 
cutting rules for loop integrals involving propagators with complex masses to show that uni- 
tarity is satisfied perturbatively jlBf. In agreement with Ref. ]2~6j . the S'-matrix connecting 
stable states only is unitary. 

We organize our perturbative calculation by applying the standard power counting of 
Refs. @, M to the renormalized diagrams, i.e., an interaction vertex obtained from an 



0(q n ) Lagrangian counts as order q n , a pion propagator as order q~ 2 , a nucleon propagator 
as order and the integration of a loop as order g 4 . In addition, we assign the order q^ 1 
to the A propagator and to the Roper propagator. Within the CMS, such a power counting 
is respected by the renormalized loop diagrams in the range of energies close to the Roper 
mass. In practice, we implement this scheme by subtracting the loop diagrams at complex 
" on- mass-shell" points in the chiral limit. 

When calculating an observable, we do not perform an expansion in powers of the mass 
differences between the Roper and the nucleon or the Roper and the A. Rather we calculate 
the chiral corrections to the magnetic moment of the Roper as a series in powers of the pion 
mass which is either divided by large scales, like AirF and the heavy masses, or multiplied 
by coupling constants which contain (inverse powers of) hidden large scales. As the omitted 
neighboring resonances, like iV(1535), couple weakly to the Roper resonance, inverse powers 
of small scales (mass differences between the Roper and the omitted resonances) which are 
hidden in low-energy coupling constants of our effective theory are enhanced by inverse 
powers of small couplings (corresponding to the weak coupling of the Roper resonance to its 
neighbors) and therefore effectively appear as large scales. 

The dressed propagator of the Roper can be written as 

i3 M = u v u\ > ( 12 ) 

f ~Z X - Z R (p) 

where — iTi R {j>) denotes the sum of one-particle-irreducible diagrams contributing to the 
Roper two-point function. The pole of the dressed propagator Sr is obtained by solving the 
equation 

z-z x -Z R (z) = 0. (13) 

We define the pole mass and the width as the real part and (—2) times the imaginary part 
of the pole [29] , respectively, 

z = m R - i — . (14) 
Some of the phenomenological analyses and dynamical models describe the Roper as a 



double-pole structure (see, e.g., Refs. |3jJ,|32j). As the self energy in Eq. ([TBI is a multi- valued 
function, one might be tempted to look for several solutions of this equation. Although the 
numbering of sheets is a matter of convention, it is our understanding that in the standard 
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FIG. 1: One-loop self-energy diagrams of the Roper. Dashed and solid lines refer to the pion and 
nucleon, respectively, and double-solid lines correspond to the Roper and delta. The numbers in 
the vertices indicate the chiral order. 




(T1) (T2) (T3) 



FIG. 2: Tree diagrams contributing to the elastic electromagnetic form factors of the Roper 
resonance. Double-solid and wiggly lines correspond to the Roper and external electromagnetic 
source, respectively. The numbers in the vertices indicate the chiral order. 



nomenclature only poles on the second sheet are relevant for the physical amplitude and 
should be interpreted as resonances. Within our perturbative approach, Eq. (TT3"|) has a 
unique solution on the second sheet. This solution is obtained as a power series in terms of 
the expansion parameter(s) of the perturbation theory. 

Close to the pole the Roper propagator can be parameterized as 

iS R {p) = ^— — + n.p. . (15) 
p - z 

The residue Zr (wave function renormalization constant of the Roper) is a complex-valued 



quantity and n.p. stands for the non-pole part. This is in full agreement with Ref. [30 
where we have shown that physical quantities characterizing unstable particles have to be 
extracted at pole positions using complex-valued wave function renormalization constants. 
Up to 0(q 3 ), Zr is obtained by calculating the Roper self-energy diagrams shown in Fig. [Q 
We do not give its explicit expression here. 



IV. MAGNETIC MOMENT 



Using Lorentz covariance and the discrete symmetries, the most general electromagnetic 
vertex of a spin- 1/2 field may be parameterized in terms of 12 Dirac structures multiplied 
by form functions depending on three scalar variables, e.g., p 2 , p 2 , and q 2 , where q — Pf —p% 
33H35J. For charged fields, the Ward-Takahashi [36j,l37| identity provides certain constraints 



among the form functions. For a stable particle such as the nucleon, on-shell kinematics 
corresponds to pf = p 2 = mjy, and the form functions reduce to conventional form factors 
of q 2 , say, Dirac and Pauli form factors F\ and F 2 , respectively. For unstable particles such 
as the Roper resonance, the analogous kinematical point is given by the pole position, i.e., 
p 2 = p 2 = z 2 . In Ref. [30| we described a method how to extract from the general vertex only 
those pieces which survive at the pole. To that end, we introduced "Dirac spinors" w i and 
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FIG. 3: Loop diagrams contributing to the elastic electromagnetic form factors of the Roper res- 
onance. Dashed, wiggly, and solid lines correspond to pion, nucleon, and external electromagnetic 
source, respectively; double-solid lines correspond to the Roper and A. The numbers in the vertices 
indicate the chiral order. 

wi with complex masses z which essentially correspond to half of the projection operators 
A + = ^jW^w^ used in Refs. 34|, |35| for the initial and final lines. In terms of these "Dirac 
spinors," the renormalized vertex function for p 2 = pf = z 2 may be written in terms of two 
form factors, 



1 Z R w t (p f )T tl (p s ,pi)w :i (p i 



w l (p f ) 



+ 



2m N 



w°(pi) 



(16) 



where is the physical mass of the nucleon. 3 Both electromagnetic form factors of the 
Roper are complex-valued functions even for q 2 < because of the resonance character of 
the Roper. As in the case of an on-shell nucleon, the third form function vanishes at the 
pole because of current conservation or time-reversal invariance. 



Note the different normalization of the magnetic form factor. 
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To 0(q 3 ), the vertex function r M (p/,p.j) obtains contributions from three tree diagrams 
(see Fig. [2]) and fourteen loop diagrams (see Fig. [3]). By multiplying the tree-order contri- 
bution with the wave function renormalization constant, one subtracts all power-counting- 
violating contributions of loop diagrams to the Fi form factor. We obtain -F\(0) = {1-\-t$)/2 
in agreement with the Ward identity. This means that, as expected, the electric charge of the 
Roper does not receive any strong corrections. On the other hand, the loop contributions to 
the magnetic form factor contain power-counting-violating terms. These parts are analytic 
in the squared pion mass and momenta. They are subtracted from the loop diagrams and 
absorbed in the renormalization of the couplings c\ and Cj. 

The anomalous magnetic moment in units of the nuclear magneton is defined as 

Kr = F 2 (0). (17) 

Since both the four-momentum as well as the polarization vector e M count as 0(q) our 
calculation yields the magnetic moment to 0(q). The tree-order result for k r is given by 

/ r* 

.tree o ™ 7 



KT = 2m N ^ + T 3 c;). (18) 

In order to show that the subtracted loop contributions satisfy the power counting we divide 
the diagrams of Fig. [3] into three separate classes. Diagrams potentially violating the power 
counting are loop diagrams with internal Roper, nucleon, and delta lines which we refer to 
as classes A, B, and C, respectively. We denote the respective contributions to the magnetic 
moment by k r B ' C . 

At first, we consider the contribution of Kr. Dividing the expression by M and taking 
the limit M ->• yields 

lim — - = -f„ tv 19 

m^oM yR 8nF 2 3 V ; 

Replacing the low-energy constant g R with g^, this expression coincides with the non- 
analytic contribution to the anomalous magnetic moment of the nucleon 0, 38 1. Next, 
we analyze the contributions stemming from k 3 ,. For a fixed and finite mass difference 
z x — m x , the limit M — > is zero 

lim -4 = 0. (20) 

M^O M 



If z x — m x scales as aM the limit M — > is given by 

k% _ 2 m N 



}}^ = 9 2 n R ^9(*)t 3 , (21) 



with 

g (a) = iir (Va 2 -I - a) +aln(2a) - Va 2 - 1 ln(a + Va 2 - 1). (22) 
Taking the limit M — > after the limit z x — > m x results in 

lim ( i im !£]=-g% R ^r 3 . (23) 

Similar results are obtained for k r . For fixed and finite mass difference z x — za x the limit 
M -»■ yields 

K C 

lim -4 = 0. (24) 

M-*0 M 
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FIG. 4: One- loop contributions to the anomalous magnetic moment of the Roper as functions of 
the pion mass. The left figure corresponds to the neutral and the right to the charged resonances. 
The solid and dashed lines indicate the real and imaginary parts, respectively. 



If z x — z Ax scales as f3M the limit M — > is given by 



lim ^ 

M 



2 m N i n\ 



Taking the limit M — > after z x — > z Ax one finds 



lim 



M->0 \ mA- 



lim -4 



M 



-9ar 



18ttF 2 



t 3 - 



(25) 



(26) 



The above analysis shows that the renormalized loop diagrams satisfy the power counting 
regardless of how the various mass differences are treated. 

To estimate the loop contributions to the anomalous magnetic moment of the Roper 
we substitute |f = 0.092 GeV, M = 0.140 GeV, m x = 0.940 GeV, z Ax = (1.210 - 
0.100 i/2) GeV, z x = (1.365 - 0.190 i/2) GeV, p = 1 GeV, g R = 1, g AR = 1, g NR = 0.45 
23[ and obtain 

KR = (0.055 + 0.090 i) - (0.223 + 0.156 i)r 3 . (27) 

Figure H] shows the loop contribution to the anomalous magnetic moment of the Roper as 
a function of the lowest-order pion mass M, where M 2 = 2Bm 0. 



V. SUMMARY 

To summarize, we have calculated the magnetic moment of the Roper resonance up 
to and including order q 3 using effective-field-theory techniques. To obtain a systematic 
power counting for energies around the mass of the Roper, we applied the CMS which is a 
generalization of the on-mass-shell renormalization for unstable particles. Unrenormalized 
contributions of loop diagrams to the magnetic moment contain power-counting-violating 
terms. However, these terms are analytic in the squared pion mass and the momenta and 
can be systematically absorbed in the renormalization of the available low-energy coupling 
constants. The renormalized loop diagrams satisfy the power counting regardless of how the 
Roper and nucleon as well as the Roper and delta mass differences are treated. 
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At next-to-next-to-leading order, 0(q 3 ), only the isovector anomalous magnetic moment 
receives a loop contribution. Analogously to the nucleon, the loop contribution to the 
isoscalar anomalous magnetic moment starts with order (9(g 4 ). 4 Due to the unstable char- 
acter of the Roper, the loop contributions to the anomalous magnetic moment feature an 
imaginary part which is of the same order of magnitude as the corresponding real part. 

At present, an extraction of the elastic electromagnetic form factors of the Roper from 
experimental measurements appears to be unrealistic. However, our expressions for the 
anomalous magnetic moment may be used in the context of lattice extrapolations. Moreover, 
lattice QCD provides for an opportunity to determine the five unkown parameters. A fit 
of our expressions to lattice data at different values for the pion mass results in a complete 
theoretical prediction of the anomalous magnetic moments of the Roper. 
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VI. APPENDIX 

Making use of dimensional repolarization with n the number of space-time dimensions, 
the loop functions are given as [42] 



B {p 2 ,m\ 



A (m>) = {2 ^t n / „ — 2 = -32ttW -2m 2 In ™ 
v ' in 2 J k 2 — m 2 + ie fi 

(27r/i) 4 ~ n r d n k 



2 2 

,m 2j - — 



in 2 J [k 2 — in 2 + ie] [(p + k) 2 — m\ + ie] 

-32tt 2 A + 21n — - 1 - - 2 F X (1, 2; 3; u) 
m 2 2 



2 \ m 2 (u — 1) J \ m 2 (uj — 1) 



u = ^ ^ VV ,\ — , (28) 



m\ — 1712 + P + V ( m i — m 2 + P ) ~ ^ m lP 

2ml 

where 2 Fi (a, b; c; z) is the standard hypergeometric function, \i is the scale parameter of the 
dimensional regularization and 



A = 

16 



^{^4i" i(te) + r ' (i > +i i}- < 29) 



4 In manifestly Lorentz- invariant baryon chiral perturbation theory, a calculation at 0(q n ), in general, not 
only produces contributions of 0(q n ) but also a string of higher-order terms of 0(q n+l ) with i = 1,2,3,... 



For the isoscalar magnetic moment, the leading term at 0(q 3 ) vanishes and only small contributions 
beyond 0(q 3 ) survive [see Eq. flU}]. 
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By writing 

F 2 (t) = m N [G 1 (t)+T 3 G 2 (t)], 
we obtain the loop contributions as: 



(30) 



Gi oop (0) 



3^ 2 



R 



- A (4)-(M 2 -34)s (4,M 2 ,4) 



\QF 2 Z X (Ad 2 - 7T 2 

V ^ H /J 64F 2 m x z x [(m x + z x ) 2 - M 2 ] n 2 

x I z x A (m 2 ) M 2 + m x (m x (m x + z x ) - M 2 ) A (m 2 ) 



M 



M 2 z 2 + m Y (m Y — z- 



in, 



+ z x y-M 2 ) b (40, 



m A 



+ (_ M 4 + m x (2m x + z x )M 2 - m x (m x - z x )(m x + z x ) 2 ) B (z\, M 2 , 



m A 



+ 



9ar 



864F 2 4 X 4 



— |-2(z Ax - z x ) (94 x - Uz x z 3 Ax + 8z 2 x z 2 Ax + 2z x z Ax + z 4 ) 



x B (z 2 , 0, 4 X ) (^a x + z x f + M 2 z 2 (94 x + ^z x z Ax - z 2 x ) M 2 

+ Az x (-64 x - 8z x 4 x + z 2 x z Ax + 4z x ) + 2 (94 x + 4z x z Ax - z 2 ) M 4 

- (^a x + z x ) (l84 x - 10z x 4 x + Tz 2 x z Ax - 34) M 2 

+ (za x + * x ) 2 (94 x - Uz x z l x + 8z l z lx + 2 4 z Ax + 4) A ° ( M2 ) 

- 2M 2 [(94 x + 4^ x 2; Ax - z 2 x ) M 2 + 2 (-94 x - 4z x 4 x + 6z 2 4 x + z\) 

x A (4 X ) - 2 [M 2 - (z Ax + z x ) 2 ] [(94 x + Az x z Ax - z 2 x ) M 4 



92 



Ax 



hz^z 



10z 2 z Ax + z x z Ax z^ 



X^Ax 



M 



+ (4 X - 4) ( 9 4 x - Uz x z l x + 8 44 x + 2z l z *x + 4)1 s ° (4' m2 ' 4 x ) f > ( 31 ) 



G l ° op (0) = 



16F 2 z x (M 2 - Az 2 x ) it 2 



-A (4) M 2 + (3M 2 - 104) Ao (M : 



+ zt 



+ 



M 2 — 2 (M 2 - Az 2 x ) Bo (4, 0, z 2 x )] - (3M 4 - lS^M 2 + 8z x ) (z 2 , M 2 , z 2 ) 



— |(3m x - % ) % A (m 2 )M 2 + 



m x (m x + z x ) 



64F 2 m x z x [(m x + z x ) 2 - M 2 } 
x (-3M 2 + 3m 2 + Az\ + 3m x ^ x ) A (M 2 ) + m x M 2 (z x - 3m x )z 2 x + m x (3m 2 + z 

M 2 - (m x + ^ x ) 2 .Bq [z 2 x , 0, m 2 ) + (m x + z x ) 3M 4 - (6m 2 + 3z x m x + 4z 2 ) M 2 



11 



r2 2 

,m x 



+ ™ x (m x + z x ) (3m 2 + z 2 x )\B (z 2 ,, M 2 

^X X ^ - 

+ 5z 2 x z Ax + 204) 



M 2 



-54zl v - 40z x z* + Q0z 2 z 2 + lOz* 



+ M 2 (27z 2 Ax + 20z x z Ax - 5z 2 x ) A (z\ x ) M 2 + 2 



54 + 20z Ax z 



+ 5 (3M 2 - hz 2 Ax ) z x - 2z Ax (lOM 2 + 534 x ) 4 _ ( 5iVf4 _ 334 X M 2 + 554 x 
+ 20z Ax (M 2 - 4 X )' + 27 (4 X - M 2 ^ Ax ) 2 A (M 2 ) - 2(z Ax + z x ) 2 
x (27z| x - 34z x 4 x - 4l44 x + 9844 x - \7z\z\ x - IOzJza* - 5z x ) 
x B (4, 0, 4 X ) - 2(M - ^ Ax - z x )[M + ^ Ax + z x ) 274 x - 34z x 4 x - 41z 2 4 x 
+ 9844 x - 1744 x - 104^ Ax - 54 + M 4 (274 x + 20z x z Ax - hz 2 ) 

- 2M 2 (27z Ax - 7z x z 3 Ax - 5044 x + 5z 3 x z Ax - 5z*\\Bo (4, M 2 , 4 X ) 1. (32) 
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